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Summary 


The Schrédinger equation is here solved with the aid of an expansion in the coupling- 
constant. Every term in the series is written as a Laplace-integral, and from this the 
asymptotical behaviour of the solution is easily obtained. The asymptotic phase will then 
also be expressed as a power series in the coupling-constant. The first term corresponds to 
the usual Born’s approximation, the higher terms can be considered as corrections for this 
formula. In this paper the explicit calculations are carried out only to the second approxi- 
mation and for small angular momenta (S- P- and D-states). The numerical results are com- 
pared with earlier approximate formulae of Pats and Ramsery and with the more accurate 
calculations made by the variational method of HuLrHEn. 


General formulae 


If the interaction potential between two particles is called U(r), the radial 
part R of the wave function satisfies the equation 


e (rR) + [e- + vin] hr 0" (1) 


For the following calculations it is more convenient to work with, instead of 
R, the function f defined by 


R=ir (2) 
and satisfying the equation 
a” P oo ee 
Es) Nea i 
(Fat) 20+ 5 + Ur 3 (3) 


As is well known, it is sufficient to consider only real solutions of this equa- 
tion, as both the real and the imaginary parts of a complex solution have to 
satisfy the equation separately. 
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We now introduce the coupling constant A, defined by 
U (r) =A (r) (4) 
where 7+g(r) is a dimensionless expression of the order of magnitude 1 (near 


the origin). 
For the solution /f we put 


f= YP fe (5) | 
n=0 


and get from the equations (5) and (3) the following formula for the com- 
putation of fy, when f,-1 is already known! 


(GE + eh) +204 1) oe + rp fn1=0. (6) 


The general solution of (6) can be written 
In = fn? + On fo (7) 


where f/f; is a particular solution of (6), Cy a constant of integration and fy 
the general solution of the homogeneous equation. The function fy is thus the 
first term in the series (5). 


Given a system of functions /;,” 


satisfying (6) we compute the function 


i Lae -> Ar foe a) (8) 


which is obviously one solution of (6). If we now use the more general ex- 
pressions (7), we get instead of f;,” the following system of functions 


n= > Cylo 
Ie me HP poe Y + Ce fo (9) 


n= - a “1 fe “ees Calg 


Equation (5) now gives the solution 
fei fy Hy Lae! fy +e = POL CA + Ca ee (10) 


. From equation (10) we see that we do not get any essentially new solution 
if we use constants Cy in (7), which are not zero. In what follows we are 
thus only interested in particular solutions of (6). 


* An expansion of this kind has earlier been used by Hrtxtunp [11]. His series is, how- 
ever, not convergent, compare FROBERG [12]. 
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Let us suppose that we know the functions fy, /,,..., fn—1 and that for the 
product r@fn—1 we-have the integral representation 


rp ina = f "hn (p) dp (11) 
Cc 


where C is some complex path of integration. 
We put for the solution of (6) 


jn = | hn—1(p) V (p) { et” W (q) dadp (12) 
Cc 


and get for the determination of W(q) and V (p) the following equation 


vo) [er + k*) W(q) + fa l—-Fe 24K) + 21a | dq=—eé". (13) 


This technique for solving the Schrédinger equation has previously been used 
by JoHansson [/]?, who has applied it to a radiation problem, and by SExt [2], 
who has solved a special scattering problem in this way. Also in the mathe- 
matical literature e.g. WHITTAKER-WATSON [3], similar methods are frequently 
used. 

We solve equation (13) by putting 


dw 


2 4 2 


+ 21qgW = P(q) (14) 


where P(q) is an arbitrary polynomial in g. Equation (14) can be integrated 
immediately, and we get 


2 
Wg) =(¢ + aya log (1 “+ 4 + B arctg 1 + R(q) (15) 


where A and B are constants and R(q) is another polynomial in g. We now 
have 


2 
a(t) = [tna (p) Vn (p) [ 14. log (1 + a + By arcte 
‘ q 


Pon ae 2\ p (16) 
er (gq? + k*)'dqdp 


2 JoHaNsson has discussed the case of a hydrogen atom with its electron originally in 
the ground state and calculated the absorbtion coefficient for radiation due to scattering. 
The problem is first treated non-relativistically, though some relativistic calculations are also 
made. It is, however, not necessary to limit oneself to the ground state, and the calculations 
of JOHANSSON can without difficulty be carried out (non-relativistically) for an arbitrary state 
of the electron. In this way one gets, after a comparison with the ordinary perturbation 
theory, a comparatively simple deduction of the polar matrix-elements of the hydrogen 
atom — even for the continuous spectrum. The results agree with those of GoRDON [Ann. d. 
Physik 2, 1031 (1929)]. 
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n+ 0, and 


fale) =o f or + #08 (1+ % Jap 
(16 a) 


z 
EOP: 


+ By | e" (p* + kh?) arctg 

The integrands in (16) and (16a) have singular points for p= +7ik. As 

we wish to obtain a real solution, we have two possible paths of integration, 
shown in Fvgs. 1 and 2. 


Hiormle Fig. 2. 
Obviously 


. 2 , 
[ Foy toe (1+ %)av = 24 | () arete® dp 


and 


2 : p 
[ (0) oe (1 =f dp =2i | / (p) aretg 7 dp 


(f{(p) supposed regular in p= + 7k). We may thus limit ourselves to con- 
sidering only the path in Fig. 2. When this contour has as passed by the 


point p, log {1 + *,) has resumed its original value and arctg 2 i has increased 


by 2a. In the first term of equation (16a) the path of integration can thus 
be modified as shown in Fig. 3. This solution is hence finite for all values 
of r. In the second term the contour cannot be simplified in this way and 
the integral diverges for r= 0. Accordingly we put By = 0 and get 
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2 
f(r) = Ay { o (p? + k*)! log (1 + 4 dp. (17) 


For large values of » we get by a well-known method (compare e.g. [3] 
p. 343). 
La 


fo(r) > —42 Aol! (2 k)'r--} sin (ir _ =) ; (17 a) 


The numerical value of the factor Ag is fixed by the condition that the ampli- 
tude is to tend to unity for large 7, and we get 
(2k)-! 


ae (17 b) 


Fig. 3. Fig. 4. 


If n is not zero in equation (16), we choose the path Cp in Fig. 4 and get 
from equation (13) 


An Vn(p) = re (p? + k2)-P1, (18) 


The constant B, in (16) corresponds to C, in (7) and according to (10) we 
put B, =O for all values of n. 
Hence we get 


fall é [in (2) (p? + yal a (q? + k?)! log [: + (‘) | dqdp. (19) 
3 


Cp 


Our next problem is to calculate the function h,(p) in equation (11), starting 
from equation (19). This is in the general case rather complicated but is an 


extremely simple problem if 
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(20) 
In this case we have 
rp (In 0) = ge f Ina (0) (9? + EY fer tla + w+ RY 
é aPeey 
2 
- log [: + aa dqdp 


The equations (17), (19) and (21) give together 


oe (4 ix) Moe alae ee 


(21) 


(p3 ae k*) 1+1 k2 


271 2 


Co.— fag 


os fee (pn + k*) log 3 oF ‘ dpn...dpy apo. 
C 


Pyn—-1 


The paths of integration are illustrated in Fig. 5. 


Fig. 5. 


It is of course possible to choose the paths OC, . .., Crn—9-u mM such a way that 


they surround more than two singular points but, as a rule, solutions ob- 
tained in this way will not be real. 
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The first (Born’s) approximation 


If n=1, we get from formula (22) 


ae [(p + w)? + Ry}! (p + p)? 
h=gz4o | Pe ra oe |e ep Ie | 
CG 


( 
i (23) 
few + k*)' log (1 + {2) dado. 
Cp 
For large values of 7 we get from this expression 
—Uttk 
cos (i = 2) 
Da ee [e+ w+ PY mn 
f(r)  S; pti (p? + yA (24) 
—U—-tk 


If in the last integral we make the following transformation of the integration 
variable 


A eG! 
2 
Pea 2ak 2) 
2 
nes ths 
B= 1 a je 
we get 
—uttik 2 ae =p | 
=f: ar a : 
i o Ar dp = 3 (—, 5} es mae 1 Q(— B) (26) 
—U—tk nl 


where Q;(—f) is a Legendre function of the second kind ([3] p. 316). Com- 
bining the equations (24) and (26) we get the asymptotic form of the func- 


tion f, (r) 


lx 
1 a cos | kr — 9 
foe da(e 2) 2 


This formula has earlier been proved by Huttuin [4] in a slightly differ- 
ent way. 


39 


G. KALLEN, The Yukawa potential 


The second approximation 


We here first discuss the S-states. In this case we get as the asymptotic 


form of (22) 


—uttk 


hope f dp 
2" 4A ker p? + ke 


=—f—tk 


To make the arctg-functions and the logarithms unique, we make a cut along 
the real positive axes of the complex p-plane and, after another integration, 


obtain 
—uUttk 
coskr 1 [#@+a2 k? ; p= dp 
fe z eal 9 log (1 + 473) + 20h [ arcty k aa 
—u—tk ' 
? = arctg 26 
lu 
The substitution e =z here gives 
—et+ik if 5 
EY cop ae hee 
Qik Ip 24 QB arctg i = 5 lon (1+ 4 | eg (og) 
: uU—2tk 
lu 
ae w+2ik 7 | 
+ ] = ] — ——_ 
| oH ( " ) ou z ian) | 
u 
u+2ik 


Putting (compare [3] p. 280) 
log (1 + & 


$ (x) = [22 De 


and 
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From the definition (32) it is easily proved that 


1 1 

6 (2) + 4 (2) = oe a} (33) 

(compare also [5] p. 344) and thus 

{ aed | ( all 
2 Re: == = lie 
oe et 2k 2 108 ee 2tk 1338) 
Noting that 
Im {¢(— ae*)} = — >) “sin (nb) (34) 
n=1 


we get after a few simple transformations of the equations (29)—(34) 


kr 1 k? = in (2nd 
pee ; (F a) lor (1+ 4%.) v2) 3 i De (35) 
r 4k 2 lu || | 
nN Laas 


= > n—* sin (nx). (35 a) 


where 


For comparatively large values of * (' aS 05) the series 


yareey o 


converges sufficiently rapidly to be useful for numerical calculations. The func- 
tion p(x) has been tabulated by CiausENn [6] to an accuracy of 16 decimals. 


If Z is small, it is more convenient to calculate the difference between (36) 


and (35a) directly, using e.g. the formula 
sin (2n0 
P= 
n? (1 5 "| 


=log K: arctg 


(37) 


K 
as sin (20) _ [. 0g ydy 
Ek coa(20) ee) 1 + y? — 2y cos (28) 
i 
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where 


ke) | 
Ea|le4 a (37 a) 


This formula can be proved if the complex integration in the equation (32) is 
done along the path shown in Fig. 6. The above difference is obviously the | 
integral along the last straight part of the path. If we here use the absolute 


value of the integration variable as an independent quantity, we get 


Fig. 6. 


K 


{2 rea VE sin (2 3) 
y = 1 — y cos (29) 


1 


An integration by parts enables us to transfer this integral to the nght-hand 
side of equation (37). If : is small (K ~ 1), the integral in (37) can easily be 


evaluated numerically. It must be observed, however, that for these small 
energies the first two terms in the series (5) give only a poor approximation 
of the wave function, unless the coupling constant is very small. 

For the asymptotic phase of the S-states we now get the following ex- 
pression 


A Le 
tg 7, = wi log (1 Se 73) a 


Vee (3 kK? VY sin (278) 
b nts Nal ee eat log Salis (G Mire eae 
AR (5 ») log (1 Sp | p (2d) + ( ay 


Beary 
we 


(38) 
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This formula can be expected to give a good approximation of the asymptotic 
phase if the energy is large and the coupling constant is small. The following 
table gives the values of tg computed with the formulae (27) and (38) com- 
pared with the values given by the variational formula of Hutruén [7]. 


k A == 0.25 0.50 0.75 1.00 1,25 1.50 1.75 2.00 
Ke “ 
Ist approx. 0.0885 0.171 0.256 0.35 0.43 0.51 0.60 0.68 
2.0 | 2nd approx. 0.0900 0.183 0.279 0.38 0.48 0.59 0.69 0.80 
numerical 0.0908 0.186 0.287 0.39 0.50 0.62 0.75 0.88 
Ist approx. 0.096 0.19 0.29 0.38 0.48 0.58 0.67 0.77 
1.5 | 2nd approx. 0.098 0.20 0.31 0.42 0.54 0.67 0.80 0.93 
numerical 0.098 0.22 0.35 0.50 0.65 0.80 0.94 1.08 
Ist approx. 0.101 0.201 0.302 0.40 0.50 0.60 0.70 0.80 
1.0 | 2nd approx. 0.105 0.219 0.342 0.47 0.62 0.77 0.92 1.09 
numerical 0.104 0.218 0.345 0.49 0.68 0.92 I Ba7l 1.85 
Ist approx. 0.0866 0.173 0.26 0.35 0.43 0.52 0.61 0.69 
0.5 | 2nd approx. 0.0934 0.202 0.32 0.46 0.61 0.78 0.96 1.16 
numerical 0.0944 0.209 0.35 0.55 0.84 133 2.40 Wels 


As a conclusion it must be said that the formula (38) gives a very poor 
approximation in those cases which are of practical interest. 
For higher angular momenta we get instead of (28) 


l 


k lz —uttik p-u p- 
i cos| = 7 [(v + | 
Je > 4 pt (p? + yA 


—u—ik Ik ik —ptik (39) 


+1 & (a:) & (a) 
oa l—7y ( f. Pi aHvy 
oa ot Seal 


where 


wot Qik | 
28 uw — Qik 


(41) 


1 ae 
ho =o Witte: yt 
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and 
2 


lu 
Se 

B 2 k2 

as before. | 
In the formula (40) the y-integrations can always be carried out with the 
aid of elementary functions, but this is not so for the z-integration. Putting © 
J; for the whole integral in (40), it can, however, be seen that the difference — 


Ji —[Pr(B)P Jo (42) 


can always be expressed in logarithms, arctg- and rational functions. (P(A) 
is here a Legendre function of the first kind). To prove this we have only to 


examine the coefficient of oe when the y-integrations have been made. 
This coefficient is 

Te ely FSi 4 1 

2 al api" a | = Py 


and hence the expression (42). 
The remaining integrations are elementary but rather tedious. For the P- 
states we get 


3 2 
J Jy = — ih-* {28 (0-7) + 4 low (1+ 5)! (43) 
and for the D-states 
2 42 
3,—-8F eae iy 
44) 
(20. bie ( ‘) Ba | ( 
. 3 1) (0 4 : og >)-—-s 
5 C= Omg) Ae OP = 4B) Naat gay 
where 
k 2 
x = arctg 7 (45) 


The second approximation of the asymptotic phases will then be given by 
es pb Ne 
tg 4p = Ok |-5 log (1 p 1] = i| - 
1. es (se 3 pee iS 
Hs | 2h eae lee (1 3 a + pr. (=<) log (1 + 17) — (46) 
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set ages RP) 3 28 


idea Il: 0) tog (1+ 4) — (28) + S aon eee 


¢ I n 
| po (1 a 1) 


Formulae useful for the computation of the asymptotic phases in the higher 
states have been given by Pais [8] and Ramsty [9]. In his paper the latter 
author has given a number of applications, and by way of comparison the same 
cases have been evaluated numerically with the use of formulae (46) and (47). 
The result is shown in the following table. 


a a 
ne - i | t2 7B tgNp | ten tg 7 pe t2 Mnum 
| | | | 
1 —5.12 | 0685 | —0.341 — 0.279 — 0.253 — 0.206 — 0.256 
1 —5.12 0.565 —0.253 | —0.209 —0.187 | —0.156 —0.190 
1 —0.84 | 0.685 —0.0560 | —0.0538 | —0.0527 | —0.0530 = 
2 2.52 | 0.685 0.0348 0.0360 0.0366 (0.0368 = 


Here tgyz and tg, are the first and second approximations used in this 
paper. Tgp, tgyr and tg num are the values given by the formulae of Pats, 
Ramsey and by numerical calculations. (These last three columns have been 
taken from the paper of RAmsry). 

Numerical calculations for the P-states have also been made by Hansson [/0] 
and we are appending a comparison with one of his results, too. 


aon gs A —0.9 2.7 43 
i Na 
Ist approx. — 0.0736 0.221 — 0.368 
2nd approx. — 0.0685 0.267 | —0.241 
Hansson — 0.0691 0.291 | —0.287 


ue, Ame 
In these applications it has turned out that the coefficient of 42 2 (46) and 


(47) is given as a difference between comparatively large numbers. In the D- 
state, for example, this coefficient is 0.000580, but comes out as the difference 
between 14.752872 and 14.752292. This means that the sums (35a) and (36) 
must be computed with a much greater accuracy than might at first be ex- 
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pected. The reason for this can be inferred if we compare the formulae (39) 
and (42). From (39) it is seen that the expansion of J; in ascending powers. 
of k begins with k?/+1, The term [P:(f)/?Jo, on the other hand, contains. 
k-2"+1, and hence the first powers of & must cancel in the end formula. © 
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